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CN ■ Abstract 

. We report results of a potential model for mesons in D space-time dimension developed 

I by considering the quark-antiquark potential of Nambu-Goto strings. With this wave 

function, we have studied Isgur-Wise function for heavy-light mesons and its derivatives 
like slope and curvature. The dimensional dependence of our results and a comparative 
study with the results of 3-1-1 dimensional QCD are also reported. 
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' ■ 1 Introduction: 

CD ■ 

^ I The two quark composite system (mesons), with one heavy quark, has been the focus of interest 
for many years [1]. In the semileptonic transitions for mesons[2], in the infinite quark mass 
^ limit, all the mesonic form factors can be expressed in terms of a single universal function, the 
(3 ! Isgur-Wise (IW) Function [3]. 

I The wave function for the heavy-light mesons has been calculated earlier within the framework 
of QCD potential model [4,5] with considerable accuracy. This has been deduced both with 
coulombic term in potential as parent and also with the linear confinement term as parent [6- 
O 9]. The characteristics like slope (charge radii ) and curvature ( convexity parameter ) of 
^ IW function has been reported in both the two cases with certain limitations. In this paper, 
we calculate the wave function of heavy-light mesons by solving D dimensional Schrodinger 
equation with potential inspired by Nambu-Goto strings. 
' Nambu-Goto action of bosonic string[10-12] predicts the quark-antiquark potential[13] to be: 



m 
in 



V{r) = ar + /io + ^ (1) 



The coefficient7 = — ^^^^j^ is the universal Liischer coefficient [14] of the Liischer term ^ , which 
depends upon the space-time dimension d. a is the string tension whose value is 0.178 GeV"^, 
fio is a regularisation dependent constant. 

With this potential, considering Liischer term as parent, we employ higher dimensional Schrodinger 
equation, for getting the unperturbed wave function. Next, we consider linear confinement term 
as perturbation and find the total wave function for mesons using Dalgarno's method of per- 
turbation. This we take as a generalisation of earlier work [15] in 3-dimension, mathematically 
considering Liischer term in Nambu-Goto potential to be identical with coulombic term in QCD 
potential. Although linear term in this potential is the leading one,we justify our such consid- 
eration under the observation that the Luscher term becomes more and more significant as D 
goes on increasing. 
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With the wave function generated, we then deduce Isgur-Wise function (IWF) and its deriva- 
tives (slope and curvature). As stated above, the same study has been carried out in three 
dimensional Euclidean space with Cornell potential [15]. In our work, we also make a compara- 
tive study of our results in higher dimension with those obtained with standard QCD potential. 
The section 2 contains the essential formahsm, section 3 , the results and the section 4 contains 
conclusion and remarks. 



2 Formalism: 

2.1 Potential Model: 

For simplification of our work and for better comparison with Cornell potential, we put the 
Nambu-Goto potential ( equation (1) ) as : 

y(r) = -^ + (7r + //o (2) 
r 



now with 7 
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Here d is the space-time dimension with d = D + 1 , D being the space dimension. 

We take — ^ as parent and ar + jiQ as perturbation. Our unperturbed Hamiltonian [16] is 
then, 

with 

H'^ar + fio (4) 
as perturbation, fx is the reduced mass of the meson. 

2.2 Wave function with Liischer term as parent: 

In D-dimension , the Schrodinger equation is [17-21]: 

[-^V^ + V^)(r)]*(r, Qn) = ^*(r, Qd) (5) 

with , 

^2 ^ J_d_,n-id_. _ ^U^d) D-l± _ Ki^n) .g) 

^ r^-^ dr dr dr^ r dr 

and 

'^{r,QD)^R{T)Y{nD) (7) 

Here, ^^^^"'^ is a generahsation of the centrifugal barrier [22] for the case of D-dimensional 
space-time. 

The eigen values of K\,{Q,d) are given by : 

kl{SlD)Y{SlD) = l{l + D-2)Y(SlD) (8) 
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Here Y^Qo) and R{r) are the spherical harmonics and radial wave function; / is the angular 
momentum quantum number and E is the energy eigen value for the unperturbed wave function. 
This gives the equation (4) in term of radial part as : 

D-ld l(l + D-2) 2a, ^ x 
For I — and taking h — 1, equation (9) yields, 



R"(r) + - — -R'(r) + 2i^(E + -)R(r) = (10) 
r r 

Following reference [23] , we look for the solution of eqn (10) , by taking the radial wave function 
of the form : 

R{r) = F{r)e-^^'' (11) 

Here, it is to be noted that the negative sign in the exponent ensures the square integrability 

at the origin and infinity. 

Using equation (11) in (10), we obtain, 

F"(r) + (^ - 2/x7)F'(r) + (^^f - + 2fiE + '^)F{r) = (12) 

Now, we take, 

00 

F(r)- J]a„rV(r,L') (13) 

n=0 

such that /(r, D) = 1 at D — 3. We take [ ] /(r, D) — , which satisfies this condition. We 
find the radial solution as: 

00 

R{r) = ^ a„r"+^e-''^'^ (14) 

n=0 

With this radial wave function we construct the wave function ^"(r) for unperturbed Hamil- 
tonian, with n = as : 

^O(r) = Nr^e-"^'' (15) 
N is the normalisation constant, which we obtain from the condition: 

poo 

/ 47rr2 I *°(r) dr ^ 1 (16) 
Jo 

Applying equation (15)in (16), we get, 

_(2^ /2 (17) 

At D = 3, this gives = (^^^)^/^ which is similar to the value obtained earlier [15] for the 
corresponding case with QCD potential. With this unperturbed wave function, we measure the 
eigen energy E as : 

E^W^^ - r ^47rr2 | *°(r) \^ dr = --f^ (18) 
Jo 1" D — 1 
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2.3 Wave function with linear term as perturbation: 

The first order perturbed eigen function ^'{r) can be calculated using the relation [22] : 

ifo^'(r) + H'-^%r) = W^°^'(r) + W^'^°(r) (19) 

where, 

poo 7~) 

' = / Anr^H' \ *°(r) \'' dr = ^ + fio (20) 
Jo 2/^7 



W 

Then, from equation (19) we get. 



Taking h = 1 and expanding, 

+ + ^ - 2^W^°]^'(r) = 2;.(ar + /.o - l^')*°(r) (22) 

Prom equation (22), following Dalgarno's method [25-27] of perturbation, we get [Appendix- A] : 

■q,'(r) = -^rV'^e-''^'- (23) 

With this perturbed wave function, we construct the total wave function as : 

^*°*«'(r) = Ni[^^{r) + *'(r)] (24) 
Using equations (15) and (23)in (24): 

^*°*«'(r) ^Ni[l- ^ry^e-"''' (25) 

where Ni is the normalisation constant for the total wave function and is obtained from : 

/"OO 

/ 47rr^ I **''*"^(r) |2 c/r = 1 (26) 
Jo 

as ^ 

^1 ^ n nz\ r(D) 2aD rp+2) a2j^2 r(g+4) n/2 (^''^ 

For D — 3, Ni gives the value of eqn. (12) of ref [15] when 7 is replaced by ^ and cr by b. 
Now, if there is only coulombic term as parent, with no perturbation, then cr — > . In this case 
eqn. (27) becomes: 

'^^-2V^[^]V2-[4^] 

This is our normalisation constant N for unperturbed wave function ^'°(r) as obtained in equa- 
tion (17). 
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2.4 Isgur-Wise Function: 

Isgur-Wisc function, ^{v,v') depends only upon the four velocities and Vi,r of heavy particle 
before and after decay. This ^{v,v') is normalized at zero recoil [29]. If y = Vi^.Vi,' , then, for 
zero recoil (y=l), ^{y) = 1. In explicit form IW function can be expressed as : 

e(y) = l-p'(l/-l) + C(y-lf (29) 



is the slope parameter at y—1 given by 



P --^\y=^ (30) 

p is known as the charge radius. 

C is the convexity parameter given by - 

The calculation of this IW function is non-perturbative in principle and is performed for different 
phenomenological wave functions of mesons [28] . This function depends upon the meson wave 
function and some kinematic factor, as given below : 



^=^1.-1 (31) 



oo 

2iiTf 1 2 



C{y)= / 47rr^|^'(r)p COS (pr)(ir (32) 







where cos(pr) = 1 — ^-j- + ^ — h with p'^ = 2p?'{y — 1). Taking cos(pr) up to 0{r^) we 

get, 

/•OO noo 2 /"CxD 

i{y)= A7rr^\^{r)\'^dr-[Anp^ /|^(r)|^(ir](2/-l) + [-7r//M r^\^{r)\'^dr]{y-lf (33) 
Jo Jo 3 Jo 

Equations (28) and (32) give us : 



p' = [4V/ r^|*(r)|'rfr] (34) 
Jo 

2 POD 

C=[-V/ r6|*(r)|'dr] (35) 

3 Jo 

47rr2|^(r)|2(ir = 1 (36) 







Equation (35) gives the normahzation constants N and N' for ^°(r) and ^*°*«'(r) as obtained 
earher in equations (17) and (27). 

2.5 Derivatives of Isgur-Wise function: 

Taking the unperturbed wave function ^°(r) from equation (17), we have calculated the slope 
and curvature of IW function ^ (y) as given below : 



2 DjD + l) 



D(D + l)(D + 2)(i. + 3) 

967^ ^ ' 
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Thus, for D = 3, wc get = and C = ^ , which are the expressions (25) and (26) of ref. 
[15], if 7 and a are replaced by ^ and b of QCD potential as given below in eqn. (38). 

An 

V{r) = -— + br + c (39) 
3r 

Also, With total wave function (r) we have calculated the slope and curvature. 



r r(D+2) _ 2aD r(g+4) , a2^2 r(D+6) i 

^ ^ r r(g) _ 2aD r(J)+2) , ^2^2 r(D+4) i 

K2/i7)0 67 (2/i7)0+2 "T 30^2 (2^^)D+4j 

, r r(D+4) _ 2crD r(D+6) , JXD+8)_] 

^ K2M7)^+" 67 (2/.7)°+^ ^ 367^ (2^7)°+" ^ /.-.X 

6 r rCD) 2ffD r(i:>+2) a2£,2 r(g+4) i v^-^'' 

L(2/i7)D 67 (2/i7)D+2 -r 36^2 (2ju7)i3+4j 



These at D = 3 becomes: 



r 3.4 o- 3.4.5.6 _|_ 0-2 3.4.5.6.7.81 



2 _ ..2 42/x7)2 7(2/^7)^ ' 47^ (2m7)6 J . s 

^ ~^ M cT 3.4 I 3.4.5.6 1 V^^^^ 

l-^ 7 (2^*7)2 472 (2^17)4 J 

4 f 3.4.5.6 _ a 3.4.5.6.7.8 , 0-2 3.4.5.6.7.8.9.IO 1 
^^/^ L(2j^7)-^ 7 (2^7)'' ~^ 472 (2^7)8 J 

fi n _ g 3.4 , <t2 3.4.5.6 1 y^-^J 

L-^ 7 (2^X7)2 4t^(2|U7)4J 

and C expressions in equations (39) and (40) also give back equations (36) and (37) when 
(7 — )■ 0. Again, at D — 3, eqns (39) and (40) transform to eqns (36) and (38) of ref. [15] on 
replacing ^ for 7 and 6 for cr , if we neglect the reletivistic effect of ref [15]. 



3 Results: 




(a) C{y) vs y for diff. D with *0(r) (b) Civ) vs y for diff. D with ^f*''*«'(r) 

Figure 1: Variation of IWF with y for B meson with diff. D values 

We have studied the variation of IW function with y for different space-time dimension D for 
B meson both for unperturbed and total wave functions taking string tension a = 0.178Gey^ 
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(fig.l). The boundary condition ^{1) = 1 is found to be satisfied everywhere. 
The variation of derivatives of IW function ( and C ) with space-time dimension D for B 
meson taking unperturbed wave function ^°(r) and total wave function ^^^i^^^r) are given in 
table 1. 

We find that with increase in D value, the slope and curvature decreases. With total wave 



Table 1: Values of normalisation constant and derivatives of IW function for B meson. 



D 


N 




C 




P' 


C 


3 


0.0047 


175.04 


12766.00 


2.49 X 10-6 


817.1 


178802 


4 


0.0031 


72.93 


1861.70 


6.5 X 10-6 


262.9 


17584.2 


5 


0.00175 


48.62 


735.50 


6.95 X 10-6 


146.18 


5218.3 


9 


0.0001186 


26.80 


175.46 


1.256 X 10-6 


hAAA 


650.2 


10 


0.00005856 


25.07 


148.56 


6.94 X 10-^ 


48.1 


498.7 


26 


5.13 X 10-^° 


17.78 


60.925 


1.106 X 10-^1 


23.67 


105.9 



function, the variation of and C with space-time dimension D for B meson is shown in fig. 2. 
As D — > oo we have also obtained the asymptotic forms of and C from equations (40) and 
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(41). 
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(a) p2 vs D with (r) 



(b) C vs D with ^^^°^<'\r) 



Figure 2: Variation of and C vs D for B meson 



asym 



P' 



2 



2cr 1 _|_ 



36(^)2 (2m^)6 



1 - 



2a 1 I 1 

6(i^) ^ 36(^)2 (2^^ 



^asym — 



2a 1 



+ 



36(^)2 (2^,^)S 
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2a 



1 



+ 



(44) 
(45) 



The asymptotic values of and C are found to be 14.5865 and 35.4608 respectively, for B 
meson. 



4 Conclusion and remarks: 

From table 1 and fig. 2, we find that, with increasing D , and C values go on decreasing 
eventually reaching the asymptotic limit. We have shown our result for B meson only as a 
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representative case. Similar pattern of result is being found to follow for other B and D sector 
heavy- light mesons like D, Dg, Bg, B^. 

Although expressions for and C give back the corresponding expressions obtained with stan- 
dard QCD potential ior D = 3 (ref. [15]) by replacing 7 with ^ and a with b, but the values 
of and C is higher in the present case for D = 3 . This is due to much lower value of 7 for 
lower D , as compared to that of the corresponding term ^ in the previous work with QCD 
potential (ref. [15]) 

However, at higher D when 7 becomes more and more dominant, our and C values go on 
decreasing. This higher dimensional approach also in turn supports our consideration of treat- 
ing Luscher term as parent with confinement term as perturbation in the present formalism. 
Lastly, we conclude by making the optimistic comment that the picture will possibly improve 
at lower values of D, if we treat linear potential term as parent with Luscher term as pertur- 
bation. Development of formalism with the latter approach is under consideration. 
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A Appendix 

From equation (22), the Schrodinger equation for first order perturbation is: 

[t^ + ^^T- + ^ - 2^W^°]vI/'(r) = 2i2{ar + /io - W')^'{r) (A.l) 
dr"^ 2 dr r 

Following Dalgarno's method [25-27] of perturbation, we obtain : 

$'(r) = (ar + /io)i?(r) (A.2) 

Taking, 

R{r) = G{r)e~''^' (A.3) 
and then, expanding G{r) in terms of power series. 



oo 

D-3 



G'(r) = ^A„r"+— (A.4) 

n=0 

we ultimately obtain: 

oo 

^'(r) = {ar + fio) ^ A^r^+^e"'^^'' (A.5) 

n=0 

This upon expansion gives- 



^'(r) = [/ioAor° + (/io^i + aAoy + (/io^a + aA^y + 

ifioAs + aA^y + (/io^ + (rA^y + ]r^e-'^^^ (A.6) 

oo 

^'(r) = CnT^r^e-^^' (A.7) 



n=0 

where. 



Co = /^o^o cind (A. 8) 

Ci+i = fioAi+i + aAi (A.9) 
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with / = 0, 1,2,3, 

Using (A.7), (A.8) and (A.9) in (A.l) , we get, 



= for I 7^ 1 and 



The equation (A.7) then reduces to: 



Now, using equation (20), viz, fj,^ — W' — — equation (A. 12) becomes 

*'(r) = [/.oAo-|^r>'^e-''^'^ 
As Aq is undetermined, we set /iqAq to be zero. This results in : 

67 
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